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Abstract 



We have analyzed the ab-initio local density approximation band struc- 
ture calculations for the family of perovskite oxides, LaMOs with M=Ti-Ni 
within a parametrized nearest neighbor tight-binding model and extracted 
various interaction strengths. We study the systematics in these interac- 
tion parameters across the transition metal series and discuss the relevance 
of these in a many-body description of these oxides. The results obtained 
here compare well with estimates of these parameters obtained via analysis 
of electron spectroscopic results in conjunction with the Anderson impurity 
model. The dependence of the hopping interaction strength, t, is found to be 
approximately r~ 3 . 



Introduction 

Electronic structures of transition metal oxides have attracted a great deal of attention 
in recent years, arising from their unusual electronic and magnetic properties. Simultane- 
ous presence of strong electron-electron interaction within the transition metal 3d manifold 
and a sizeable hopping interaction strength between the 3d and oxygen 2p states is pri- 
marily responsible for the wide range of properties exhibited by these compounds. Often 
the presence of a strong intraatomic coulomb interaction strength makes a single-particle 
description of such systems inadequate, necessitating a model many-body Hamiltonian ap- 
proach mi . However, such parametrized approaches require a prior knowledge of the various 
electronic interaction strengths, such as the intraatomic Coulomb interaction strength, Udd, 
the charge-transfer energy, A, and the metal (i-oxygen p hopping interaction strength, t, 
in order to provide a description of the ground state electronic structure. Traditionally, 
various high-energy spectroscopic results in conjuction with model Hamiltonian calculations 
have provided estimates for such interaction strengths 0-§J . However, it is well-known that 
wide ranges of parameter values are compatible with experimentally obtained spectra, lead- 
ing to non- unique solutions 0. This arises from the fact that the effect of a change in 
one parameter on the calculated spectrum can often be compensated by suitably changing 
another parameter, such that the the various solutions are indistinguishable within experi- 
mental uncertainties. This problem has given rise to different sets of estimates of interaction 
parameters from different groups for the same compound. In this context, it is highly de- 
sirable to obtain independent estimates of various interaction strengths from methods other 
than those using spectroscopic results. Even if only one interaction strength can be reliably 
estimated from any other method, the problem of non-uniqueness in analysing experimen- 
tal spectroscopic result is considerably eliminated and such an approach will yield a more 
consistent description of the electronic structure. 

Recently we have shown || that the electronic and magnetic structures of a particular 
class of transition metal oxides, namely the 3d transition metal perovskite oxides of the gen- 
eral formula LaM03 (M=Cr-Ni) are described very accurately within ab-initio approaches 
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based on local density approximation (LDA). Futhermore, it is also seen that the various 
excitation spectra, such as x-ray photoemission (XP) spectra of the valence band region and 
bremsstrahlung isochromat (BI) spectra of these compounds are also described well within 
the band structure approach. This success suggests that such band structure results could 
be useful in providing reliable estimates of various interaction strengths in this interesting 
class of compounds which have attracted a lot of attention in recent times ||. Estimates 
of various interaction strengths based on ab-initio calculations have indeed been carried out 



in the past, for example for the 3d transition metal monoxide series [jn|, La2Ni04 |TT 



La2Cu04 |12[], the Ad transition metal oxide, Sr2Ru04 |T3[ and Mn doped CdTe [fL4}| . We 
report here the results of our analysis for the entire series of compounds belonging to the 
LaMC>3 family with M=Ti-Ni. Wherever possible, we compare the results obtained from 
the present approach with those obtained by analysis of various spectroscopic results already 
existent. 
Methodology 

We have performed detailed band structure calculations within the linearised muffin-tin 
orbital (LMTO) method using the atomic sphere approximation (ASA) [ftlj] . Spin restricted 
calculations were performed for the real crystal structures experimentally observed. Thus, 
no attempt has been made to determine the lattice constants from total energy calculations. 
This is justified in the present context, since our primary aim here is to obtain strengths 
of the electronic interactions in the real systems with the observed lattice constants. In 
the case of LaCo03, however, we have performed calculations for various values of the 
lattice parameter in order to deduce the dependence of the hopping interaction strengths 
on distances. Analysis of these results shows that the present method predicts the lattice 



parameters with nearly 1% accuracy. LaTiC>3 UM, LaV03 ||16||, LaCr03 |T7| and LaFe03 



I8| in the observed Pbnm structure were calculated with twenty atoms in the unit cell, 



while LaCo03 [|T^] and LaNi03 []2D] calculations were performed with ten atoms in the unit 
cell in the real R3c structure. In the case of LaMn03, the real crystal structure is Pnma 
stabilised by a Jahn- Teller distortion around the Mn-ions ||21|| . In contrast to all other 



LaMC>3 compounds, this leads to a distorted octahedra of oxygens with two distinct Mn- 
O distances. Since, this would introduce additional parameters in the analysis, we have 
preferred to calculate the band structure of LaMnOs with an earlier reported R3c idealised 
structure, where the Jahn- Teller distortion is suppressed [p^j . Such approximations have 



been made in the past j33[ and we do not expect to introduce any significant error in the 



parameters obtained, since the change in the interaction parameters is expected to be small 
for the small distortions observed. 

Various interaction strengths can be estimated by mapping the band dispersions obtained 



from a parametrized tight-binding model onto those from ab-initio band calculations |10 
Since, we are primarily interested in the states arising from transition metal d-oxygen p 
interaction which dominate the occupied density of states as well as parts of the unoccupied 
states near the Fermi energy, E^, the tight-binding Hamiltonian consists of the bare energies 
of the transition metal d (e^) and oxygen p (e p ) states and hopping interactions between all 
these states. The nearest neighbor hopping interactions were expressed in terms of four 
Slater-Koster parameters, namely ppo~, pp^-, pda and pdn [2"1||. It is further necessary [25 



to include an interaction, sda, between the transition metal d and oxygen s states in order 
to simulate the splitting between the t 2g and e g states at the r point, since the degeneracy 
of these two states is not lifted by the p — d interaction at this symmetry point. In general, 
it is also necessary to include the effect of non-orthogonality of atomic functions located at 



different sites by considering the overlap matrix in such tight-binding approaches ||25|| . On 
the other hand, the tight-binding part of any parametrised many-body Hamiltonian, such as 
the Hubbard Hamiltonian, ignores the overlap matrix with an assumption of orthogonality of 
the basis functions. Since we are eventually interested in obtaining estimates of interaction 
strengths that enter such many-body Hamiltonians as parameters, we have carried out the 
fitting of the LMTO dispersions within two separate tight-binding models, one with and the 
other without the assumption of orthonormal basis functions. 

Since we take into account only the transition metal d and oxygen p states, there are 56 
bands to be fitted in the Pbnm structure and 28 bands in the R3c structure. LMTO-ASA 



results indicate that several of the low-lying bands within the nominal oxygen p-bandwidth 
have non-negligible contributions from La derived states. Thus, we have not taken into 
account some of these bands in the fitting procedure. We included the top 38 bands in the 
case of the Pbnm structure and 19 in the R3c structure. In the case of LaTi0 3 , the top 8 
bands with mainly Ti d character overlap extensively with other bands having primarily La 
character. Thus we could not include these bands in the fitting procedure. We also checked 
the effect of including only the primarily transition metal (i-derived bands, (top 20 in Pbnm 
and 10 in R3c structures), leaving out all the primarily oxygen p-derived bands, from the 
fitting procedures. However, in the case of LaTiC>3 we have not carried out such an analysis 
due to the extensive overlap of Ti d and La derived bands mentioned above. Results of these 
different procedures, namely assuming the orthogonality or non-orthogonality of the atomic 
basis functions and employing different limited sets of bands, are consistent with each other. 
Results and discussion 

Total density of states as well as partial densities of Mn d and O p states in LaMn0 3 states 
obtained from LMTO-ASA calculations are shown in Fig. 1. We show the corresponding 
band dispersions along various symmetry directions in Fig. 2a. The zero of the energy scale 
corresponds to the Fermi energy, E F . A set of two bands can be seen close to 3 eV at the r 
point and dispersing to higher energies along both F — L and F — Z directions, becoming 
larger than 4 eV halfway to the zone boundaries. These two bands are primarily La derived 
with negligible Mn d contributions (see Fig. 1) and we do not discuss these any further 
here. Next, one can observe a group of four strongly dispersing bands in Fig. 2a between 
about 0.5 and 3.5 eV giving rise to the approximately rectangular DOS marked A in Fig. 1. 
These bands are derived primarily from Mn <i(e 9 )-oxygen jo-admixture with dominant Mn d 
character. Below the e g bands, there is a group of 6 bands spread between -1 and 0.5 eV 
with considerably less dispersion than the e g bands (Fig. 2a) giving rise to a narrow DOS 
marked B in Fig. I. Various partial densities of states in this energy range suggest that this 
group has primarily Mn d character with some finite oxygen p admixture (Fig. 1). From 
these observations it is clear that this group of bands arises mainly from Mn d (t 2g ) states. 
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Both e g and t<i g bands discussed so far arise from Mn d-oxygen p antibonding combination, 
with dominant Mn d character. Thus, these bands are termed antibonding e* and t% bands 
respectively. At still lower energies, we find an energy region with DOS between about 
-2.7 and -4.5 eV which is almost completely contributed by oxygen p-character with very 
small Mn d character. These features in the DOS are marked as C in Fig. 1 and arise 
from essentially oxygen-oxygen interactions with a non-bonding character with respect to 
Mn d- O p interactions. The corresponding dispersions of the eight bands can be observed 
in Fig. 2a. Dispersions of the remaining ten bands can be seen between -4 and -7.5 eV. 
The corresponding density of states (marked D) peaks at about -5.5 eV and has primarily 
oxygen p-character (see Fig. 1). However, finite Mn d admixture is also observable in this 
energy range. Thus, these states arise from oxygen p-Mn d bonding interactions and are 
the bonding counterparts of the Mn rf-dominated antibonding e* and t^ bands at higher 
energies. These interpretations are consistent with the description of bonding in LaMOs 



compounds which have been discussed in detail in recent times [p6|-[28|. 

In order to map the band dispersions obtained from LMTO-ASA calculations onto the 
tight-binding model, we have used all the ten primarily Mn d-dehved e* and t* 2g bands as 
well as the top nine oxygen p-bands, as explained in the previous section. The best fit 
tight-binding dispersions are shown in Fig. 2b for comparison with the ab-initio calculated 
dispersions in Fig. 2a. It is evident from Fig. 2 that good agreement is obtained between 
the LMTO and tight-binding results; this is particularly true of the bands related to the t\ g 
and e* distributed between -1 and 3.5 eV. This is significant, since the lower bands related 
primarily to oxygen p-states are completely filled in all these oxides and the electronic and 
magnetic properties in these cases are controlled entirely by the t\ g and e* bands. Moreover, 
the present result suggests that the electronic structure of these transition metal oxides can 
be well described in terms of models involving only the transition metal d and oxygen p 
states (i.e d — p models), as is the usual procedure. 

While the above mentioned case of LaMnOs was calculated within the R3c structure, 
we show an example of LMTO-ASA DOS and band dispersions for the Pbnm structure in 



Figs. 3 and 4a respectively for LaFe0 3 . Once again, four groups of DOS features are easily 
identified in Fig. 3 for LaFe0 3 . These are marked A through D and have the same origin 
as discussed in the case of LaMnOs shown in Fig. 1. In the case of the Pbnm structure of 
LaFeOs however, there are twenty atoms in the unit cell, thereby doubling the number of 
bands in comparison with the previously discussed R3c structure. Thus, within the d — p 
part of the electronic structure alone, there are fifty-six bands; eight related to the e* bands, 
twelve to the t^ g bands, sixteen related to the so-called non-bonding oxygen p-parts and 
twenty related to the bonding states of the Fe d-0 p interactions. For the sake of clarity, we 
show the dispersions of those bands which are related to the e* and bands with primarily 
Fe d character along various symmetry directions in Fig. 4a. The strongly dispersing group 
of eight bands (more easily recognisable along the r — R direction) between about 0.5 eV 
and 3 eV arise from the e* bands, while the considerably more flat twelve bands between 
-1 and 0.5 eV comprise the t* 2g region. These band dispersions alongwith the sixteen non- 
bonding oxygen p bands were fitted with d — p only tight-binding model along all symmetry 
directions in the same way as in the case of LaMn0 3 . The resulting best fit tight-binding 
results are shown in Fig. 4b. Comparison of Figs. 4a and b indicate that the d — p only 
nearest neighbor tight-binding model describes the band dispersions observed within the 
ab-initio calculation quite well. 

The parameters for the best fits to the LMTO band dispersions for all the LaM0 3 
compounds with M=Ti-Ni are given in Table I. While the main entries in this Table are 
for the cases where d and p bands were fitted, the numbers in parentheses are obtained by 
fitting only the e* and t\ g bands. These two sets of estimates are very similar exhibiting 
some systematic changes across the transition metal series. We have plotted the strengths of 
various hopping interactions across the 3d transition metal series in Fig. 5. It is reasonable 
to expect that the variations in the hopping interactionwill be related to the changes in the 
relevant atomic distances. Thus, we have also shown the experimentally observed oxygen- 
oxygen (ro-o) and the metal-oxygen {tm-o) distances in the inset of Fig. 5. Both ro-o and 
tm-o appear to be larger in LaFe0 3 compared to the overall trend. This may arise from 
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the stability of the half-filled high-spin d 5 ionic configuration of Fe 3+ in this compound. 
The decrease of these distances for LaCo0 3 , just after LaFe0 3 , is clearly related to the 
low-spin configuration of the Co 3+ ion. Fig. 5 shows that the magnitude of ppa does not 
show any significant variations between V and Mn, then decreases for Fe followed by a 
substantial increase for Co. These variations in ppa strength can be easily related to the 
changes in the nearest neighbour oxygen-oxygen distances (jo-o) in the LaM0 3 compounds; 
ro~o substantially increases for LaFe0 3 giving rise to the observed decrease in ppa in this 
compound. For LaCo0 3 , r -o is seen to be the smallest in the LaM0 3 series and for LaTi0 3 , 
ro-o is the largest; consequently the corresponding strengths of ppa are the largest and the 
smallest respectively in the series. Various transition metal d-oxygen orbital interactions, 
namely sda, pda and pdn also exhibit systematic variations across the series, the qualitative 
behaviour of these three interactions being quite similar to each other (Fig. 5). These 
interactions, besides being influenced by the transition metal-oxygen distance (jm-o) shown 
in the inset, are also influenced by the spatial extent of the transition metal d orbitals. It 
is well-known that the d-orbitals contract across the transition metal series. Thus, the 
decreasing trends in the hopping interaction strengths between the transition metal d and 
oxygen orbitals between Ti and Mn, inspite of a small decrease in r M _ G arise from the d- 
orbital contraction across the series. Interestingly, the interaction strengths are larger for 
LaCo0 3 than LaFe0 3 . This is clearly related to the pronounced decrease of r M _ between 
LaFe0 3 and LaCo0 3 , associated with the low-spin configuration in LaCo0 3 in contrast to 
the high-spin configuration in all the other LaM0 3 compounds with M=Ti-Fe. 

Besides the various hopping interaction strengths, the bare energy difference, (e^ — e p ), 
between the transition metal d and oxygen p orbitals exhibits a monotonic decrease across 
the LaM0 3 series, as shown in Fig. 6. This is indeed the expected trend, since the transi- 
tion metal d level is increasingly stablized with respect to the oxygen p-orbital due to the 
increasing nuclear potential with increasing atomic number of the transition metal ion. 

The results summarized in Table I and Figs. 5 and 6 are obtained by fitting the LMTO 
band dispersions to the results of the tight-binding model with finite overlap between the 
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orbitals at different atomic sites. However, most of the model many-body Hamiltonian 
approaches assume an orthogonal basis set; thus, the results in Table I cannot be directly 
used to estimate the parameter strengths appearing in such models. In order to provide 
estimates for such interaction strengths that parametrize the many-body calculations, we 
have also fitted the LMTO band dispersions with the results of a tight-binding Hamiltonian 
assuming an orthogonal basis. Thus, the tight-binding model corresponds to the one-electron 
part of the multiband Hubbard model involving all the transition metal d and oxygen p 
orbitals, suitable for the LaM03 series. The resulting estimates of the various interaction 
strengths corresponding to the best simulation of LMTO band dispersions are given in Table 
II. A comparison of Tables I and II show that the estimates of various parameters are quite 
similar in the two cases, justifying the assumption of an orthogonal basis in describing 
the electronic structures of these compounds. This is further supported by the fact that the 
orbital overlaps required to optimize the simulation of LMTO band dispersions are generally 
very small. We also find that the parameter values summarized in Table II exhibit similar 
systematic trends across the transition metal series as those in Table I (see Figs. 5 and 6). 

In order to verify the relevance of the interaction parameters thus estimated, we note 
that there have been several attempts in the past to estimate many of these from various 



high-energy spectroscopic results in conjunction with different many-body approaches |29 
It is well known that such estimates are often non-unique f7j; however it is believed that 
in particular, various hopping interaction strengths can be estimated with a fair degree of 
accuracy from such approaches. Thus, we compare the estimates of the transition metal d- 
oxygen p interaction, pda obtained from high-energy spectroscopy and the present approach. 
Two different groups have systematically worked on obtaining these experimental estimates 
for LaMOs compounds. In Fig. 7, we have plotted the experimentally obtained estimates 
of pda from the Tokyo group || and the Bangalore group against the estimates obtained 
in the present work. We show the results from these two groups using different symbols; we 
have drawn two broken lines as guides to the eye for the average overall variation of pda 
across the series, as obtained by each of the two groups. The solid line in Fig. 7 with a unity 



slope and no intercept, represents the behaviour of pda, if the experimental estimates and 
the present estimates were identical in every case. Obviously, the experimental estimates 
are somewhat different from the calculated ones. However, it is clear that the average vari- 
ations of pda across the series obtained by these groups are very similar to that suggested 
by the present calculations; this is indicated by the fact that the broken lines represent- 
ing the average variation of experimentally obtained pda are approximately parallel to the 
solid line within the accuracy of experimental estimates. Moreover, while the estimate of 
pda from one group is higher than the calculated ones, the estimate from the other is con- 
sistently lower. This arises from the non-uniqueness of the parameter strengths estimated 
from experimental results discussed earlier. More importantly, it is obvious that pda val- 
ues estimated in the present work are also compatible with the high-energy spectroscopic 
data, since values both larger and smaller than the calculated values provide satisfactory 
description of the experimental observations. Thus, it is desirable to constrain the strength 
of pda to the values calculated here on the basis of ab-initio approaches; this will avoid 
the problem of non-uniqueness normally encountered in such analysis J7| and will help in 
obtaining considerably better defined estimates of other parameters in the model. 

The estimates of — e p obtained by analysing the LMTO results can be related in an 
approximate way to the charge transfer energy, A defined [HJ within the Anderson impu- 
rity Hamiltonian which is often used to provide a many-body description of the electronic 
structure of these oxides. For this, we first note that the eigen values, and e p in an LDA 
calculation are not directly related to the orbital ionization energies. In order to see the 
relationship between the charge transfer energy and the eigen values, we note that the total 
energy E within the LDA calculation can be expressed as a Taylor series expansion in terms 
of the electron occupancies, nj, of the various levels i: 

E{...m...) = E + ^2birii + ^a ij n i n j + ... 

i ij 

Retaining terms upto the second order in the above expansion, it is easy to show [25] that 



the coefficients related to various electron-electron interactions. The atomic orbital 
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eigen value e« of the level i is related |3(| to the total energy as 



dE _ 
drii 



i 



Then the charge transfer energy, A, defined as 



A = E(...n d + l...n p - 1...) - E(...n d ...n p ...) 



is readily given by 



A = e d - e p + -U dd 



(1) 



Thus, A and the difference in the orbital eigen values, e d — e p are linearly related with 
the difference being half of the intraatomic Coulomb strength within the 3d electrons. In 
Fig. 8, we have plotted the values of A obtained from various high energy spectroscopic 
results as a function of the orbital energy difference, e d — e p , calculated in the present work. 
The straight line with a unity slope represents the expected dependence of A on e d — e p in 
absence of Coulomb interaction. With the sole exception of the case of V, all experimental 
estimates of A is larger than that represented by the straight line; this clearly indicates 
the presence of finite Coulomb interaction strengths in all these cases. We however would 
not like to estimate U dd from the difference between the experimentally estimated A values 
and the value (represented by the straight line in Fig. 8) expected in absence of U dd , since 
the uncertainty in estimating A from experiment can be very large (as large as ± 2 eV). 
This is also possibly the reason why the A value in the case of the vanadium compound has 
been estimated lower than the straight line. The present results suggest that the estimate 
of e d — e p in this work can be used to guide the choice for A via the equation 1. Essentially, 
equation 1 gives a relationship between A and U dd , since e d — e p is estimated here, this leaves 
only one parameter (U dd or A) to be determined from experimental results. 

On many occassions, it is desirable to know how pda changes with the distance between 
the transition metal and the oxygen atom in a compound. This may be very useful in describ- 
ing the electronic structure of any compound under pressure, or in the case of compounds, 
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such as LaMnOs, where various transition metal-oxygen distances are observed. Thus, there 
have been many attempts to express the functional dependence between pda and tm-o- It 
was shown earlier that the interatomic matrix-elements are supposed to scale with distance 
as \/r l+l ' +1 [5^1, where 1 and 1' are the angular momenta of the orbitals involved. Thus, 
for p — d interactions, the matrix elements scale as 1/r 4 . In order to obtain an ab-initio 
estimate of the functional dependence, we have calculated the band structure of LaCoOs 
with different lattice parameters. The different sets of band dispersions have been analysed 
in terms of the same tight-binding Hamiltonian with overlap included. The parameters in 
each case were obtained by the previously discussed least-squared error fitting procedure. 
We have plotted the pda thus obtained as a function of r Co -o f° r the eight different calcu- 
lations in a logio-logio plot. The smallest r Co -o corresponds to a 4% contraction while the 
largest value corresponds to a 9.9% expansion over the equilibrium rc -o- Over this entire 
range, logio (jorfo") is found to vary approximately linearly with logio ijCo-o) with a slope of 
-3.04±0.02. This implies that pda varies approximately as l/r^t _ which in other terms is 
\/r l+v in contrast to the earlier expectation of \/r l+v+l dependence. 

In conclusion, the electronic structure of the LaM03 series have been studied by map- 
ping the results of a tight-binding Hamiltonian with nearest-neighbor interactions onto the 
ab-initio band structure results for M=Ti-Ni. All the interaction parameters are found 
to exhibit systematic changes across the transition metal series. The estimates of various 
interaction strengths compare well with estimates obtained from spectroscopic data, wher- 
ever available. The hopping interaction strength, pda is found to depend on the transition 
metal-oxygen separation, r approximately as r~ 3 . 
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I. FIGURE CAPTIONS 



Fig.l The LMTO DOS for LaMn0 3 (solid line) alongwith the partial DOS for the Oxygen 
p states (dotted line) as well as for the Mn d states (dashed line). 

Fig. 2 (a). The LMTO band dispersions, and (b) the best-fit TB band dispersions for 
LaMn0 3 along the symmetry directions TL, LF, FZ, ZF . 

Fig.3 The LMTO DOS for LaFe0 3 (solid line) alongwith the partial DOS for the Oxygen 
p states (dotted line) as well as for the Fe d states (dashed line). 

Fig. 4 (a). The LMTO band dispersions, and (b) the best-fit TB band dispersions for 
LaFe0 3 along the symmetry directions rX, rY, FZ and FR. 

Fig. 5 The variations of the interaction strengths sda, ppa, ppn, pda and pdw across 
the transition metals series, Ti-Ni. The inset shows the interatomic distances between the 
metal-oxygen (r M _ ) and the oxygen-oxygen {r -o) sites in Afor the series. 

Fig. 6 The variation of e p -e d across the transition metals series, Ti-Ni. 

Fig. 7 A comparison of the calculated (pda) ca i c to the values estimated from experiments, 
(pda) expt by different groups (Tokyo and Bangalore). The solid line has been drawn for 
(pdcr) ca i c equal to (pda) expt and the dashed lines indicate the trends seen in the data from 
different groups. 

Fig. 8 A comparison of the calculated value of e d — e p with the experimentally obtained 
A. The solid line represents A=e^ — e p . 

Fig. 9 The variation of logi (pda) with \og w (rc -o) in LaCo0 3 . 



16 



TABLES 



TABLE I. Estimates of tight-binding parameters by a least-squared-error fitting of LMTO 
results including orbital overlaps. The numbers in parentheses have been obtained by fitting 
only the transition metal d bands. 
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TABLE II. Estimates of tight-binding parameters by least-squared-error fitting of LMTO 
results with the assumption of orthonormal basis functions. The numbers in parentheses 
have been obtained by fitting only the transition metal d bands. 



Compound 
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Fig.1 
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